Cell Reports Methods

Cryo-electron tomography (cryo-ET) enables near-native visualization of subcellular architectures, yet applying it to moderately thick, multilayered tissues such as the retina is hampered by inadequate vitrification and inaccurate depth-targeting. Here, we developed PLCT, an integrated approach combining modified high-pressure freezing, cryo-ultramicrotome trimming, and plasma-based cryo-FIB milling to overcome these barriers. PLCT reliably vitrified <100 m retinal strips with minimal ice artifacts, navigates precisely to the outer plexiform layer using morphological landmarks, and produces high-quality lamellae suitable for high-resolution cryo-ET. Subtomogram averaging (STA) analysis identified microtubules at 16.33 [A] within retinal horizontal cell processes. Importantly, STA also resolved a 10-nm-diameter filamentous structure at 24.81 [A] in the same processes, featuring six peripheral strands surrounding an elongated central density with continuous intervening cavities, an architecture consistent with intermediate filaments. Together with its native localization and immunoreactivity, these features collectively identify the filaments as neurofilaments. Separately, 3D reconstruction of synaptic ribbons uncovered a previously unrecognized "mahjong tile"-like fine ultrastructure. These results demonstrate that PLCT-produced lamellae are of sufficient quality to support structural analysis in native tissue. Although demonstrated on retinal photoreceptor synapses as a proof-of-principle, PLCT is inherently generalizable, with its depth-navigation and vitrification strategies directly applicable to any multilayered tissues. This work establishes PLCT as a robust, reproducible platform for depth-resolved in situ cryo-ET of multilayered tissues.

Current approaches for imaging fibrotic remodeling have sensitivity, specificity and cost drawbacks that limit both preclinical research and clinical diagnosis. Here, we show that fast green FCF, a small molecule that binds to fibrillar collagen, enables highly sensitive and specific imaging of fibrosis in lung samples from mice and humans using fluorescence microscopy. We report strategies for using fast green FCF staining to assess fibrotic remodeling using precision-cut lung slice and whole-biopsy preparations. Our findings demonstrate that fluorescence imaging of fast green FCF-stained collagen will be useful for fibrosis research and may help to improve detection of fibrosis in clinical pathology.

Motivation: Cell-free DNA methylation provides a minimally invasive signal for early cancer detection and tissue-of-origin prediction. Most methods represent methylation measurements as independent fixed-window features and therefore do not explicitly model relationships among genomic regions. Results: We developed PANGEM (Pan-cancer Graph-based Cancer Detection Using the Cell-free DNA Methylome), a graph-learning framework that represents genomic bins as nodes and integrates CpG context, genomic proximity, and sample-specific methylation similarity in the graph topology. Across five repeated stratified train-test splits, PANGEM achieved the highest mean performance among evaluated methods, with an AUROC/AUPR of 0.997/1.000 for binary cancer detection and macro-AUROC/AUPR of 0.977/0.870 for multiclass tissue-of-origin prediction. In the independent INSPIRE cohort, 72 of 78 cancer cases (92.3%) exceeded the binary classification threshold, and PANGEM correctly classified 9 of 17 head and neck cancer cases (52.9%), the highest accuracy among evaluated methods. Subnetwork analysis further identified recurrent, graph-connected methylation patterns, including a 111-DMR subnetwork with increased methylation in cancer samples.

Meso2EM is a correlative light and electron microscopy workflow that transfers neurons selected from mesoscale functional images to targeted electron microscopy. We recorded Ca{superscript 2} signals from layer 2/3 neurons across a contiguous 3 x 3 mm cortical field in awake mice and reidentified a selected neuron after fixation and tangential sectioning. Lectin-labeled vascular architecture served as a shared landmark across in vivo two-photon imaging, confocal microscopy, laboratory micro-CT of resin-embedded tissue, and block-surface scanning electron microscopy, guiding focused-ion-beam scanning electron microscopy to the target cell body. The same progressive-targeting principle also supported serial ATUM-SEM reconstruction of an in vivo-tracked dendrite and serial transmission electron microscopy of optically selected dendrites from a patch-clamp-recorded Martinotti cell. Meso2EM therefore provides a practical route for preserving target identity across large changes in scale and specimen state while restricting electron-microscopy acquisition to a selected region.

Astrocyte border formation is a conserved neuroprotective response to neural tissue disruption, yet astrocyte border states at implanted biomaterials remain less well characterized than injury responses. Here, we developed the Astrocyte Border Characterization (ABC) Tool, which leverages a shear-thinning, injectable biomaterial to locally deliver astrocyte-specific RiboTag AAVs and small molecule regulators in the mouse striatum, enabling molecular profiling and phenotypic modulation of astrocyte border (AB) cells. Spatially precise delivery of AAV using the ABC Tool yielded enhanced specificity and robust RiboTag expression in AB cells from 7-70 days post injection. Temporal transcriptomic profiling of AB cells revealed predominantly acute, transient changes in genes governing dedifferentiation, proliferation, metabolic reprogramming, and inflammation regulation. Persistent changes accounted for only 14% of regulated genes but involved critical gain of functions in immune regulation and host defense that mirrored astrocyte border responses at chronic CNS injuries. Local delivery of indiscriminate or astrocyte-selective ablation molecules delayed, rather than prevented, border formation, ultimately yielding thicker astrocytes borders with increased inflammation and fibrosis at the biomaterial-tissue interface. Conversely, local delivery of {beta}-hydroxybutyrate (BHB) from the ABC Tool altered key aspects of the transcriptional reprogramming to attenuate chronic astrocyte reactivity and prevent biomaterial contraction without exacerbating inflammation or fibrosis. Our findings establish the ABC Tool as a bioassay for studying and manipulating astrocyte borders at implanted biomaterials and identify focal metabolic regulation as a strategy to modulate AB cell phenotypes and enhance the CNS biocompatibility of biomaterials.

Accurate, reproducible interpretation of kidney allograft biopsies is critical for diagnosis of graft injury to guide prognosis and management. The international Banff classification is a consensus diagnostic system based on semiquantitative histological lesion scoring on either extent or severity of kidney transplant biopsies. However, pathologist scoring is limited by substantial interobserver variability, constrained scalability, and the inherent nature of the scoring system itself. Here we present BanffNET, a weakly supervised, probabilistic deep learning framework that combines self-supervised feature extraction with a novel Bayesian multiple-instance learning framework to predict (continuously) the full spectrum of Banff lesion scores directly from whole-slide images (WSIs). Using lesion-specific aggregation functions tailored to localized (modeling lesion severity) and diffuse pathologies (modeling lesion extent), BanffNET generates interpretable, patch-level probability maps and calibrated slide-level scores. BanffNET's performance was assessed relative to consensus, biological correlates of rejection and clinical outcome, demonstrating superior consistency, transportability and generalization. Trained on 7,249 WSIs from three cohorts, BanffNET demonstrates consistent performance on 11,028 WSIs across five external test sets, performing on par or exceeding expert consensus across lesions. BanffNET scores align more closely than pathologist Banff scores with molecular profiles of rejection, offering a transparent, biologically grounded framework for computational pathology with relevance beyond transplantation.

Spatial omics technologies resolve molecular expression and spatial architecture at single-cell resolution, but profiling whole slides remains costly. In practice, only a few regions of interest (ROIs) are profiled, leaving the rest of the tissue unmeasured. S2-omics was the first framework to unify ROI selection with out-of-ROI prediction, but it operates on superpixels rather than individual cells and predicts discrete cell types rather than continuous molecular profiles. Superpixel-based representations do not explicitly preserve cell boundaries, while categorical cell-type labels cannot quantify molecular expression within cells. Here we present RECON, a two-stage framework that performs ROI inference and whole-slide molecular reconstruction at single-cell resolution, predicting both continuous molecular profiles and discrete cell-type labels. In the first stage, RECON extracts morphological and microenvironmental features from individual cells to identify a representative ROI for spatially resolved single-cell molecular profiling. In the second stage, RECON trains deep learning models on molecular measurements acquired within the selected ROI and reconstructs transcriptomic or proteomic profiles for all remaining cells on the slide. Benchmarked against pathologist annotations, RECONs ROI selection outperforms the superpixel-based S2-omics approaches (IoU: 0.75 versus 0.64). For transcriptomics, refining the modeling unit from superpixels to single cells improves per-gene Pearson correlation by 22%. For proteomics, RECON surpasses the current state-of-the-art method, ROSIE, across all 16 markers, with a median per-cell Pearson correlation of 0.91 versus 0.84. Moreover, RECON delineates tumour boundaries and regions with distinct immune-cell densities, and highlights candidate tertiary lymphoid structures. Together, these results demonstrate that RECON enables informative ROI selection and whole-slide molecular reconstruction at single-cell resolution for both spatial transcriptomics and spatial proteomics.

Investigating the neural mechanisms of social group interactions and other naturalistic behaviors in small animals remains limited by current technology. Tethered neural recording systems are incompatible with many of these behaviors, while existing wireless devices for small animals are constrained by weight, bandwidth, recording duration, and the lack of closed-loop modulation capabilities. To overcome these limitations, we developed a Wireless, Interactive, Lightweight Datalogger (WILD) with integrated flexible neural probes, optogenetics, an inertial measurement unit, an ultrasonic microphone, and a head-mounted camera. This platform enables simultaneous, long-term recording of neural activity, locomotor variables, vocalizations, and eye movements from groups of freely moving mice in both laboratory and outdoor settings. Model-based real-time signal processing detects specific neural events and behavioral motifs to trigger closed-loop neural interventions. By combining multimodal recordings with advanced onboard signal-processing capabilities in a compact device, WILD enables the investigation of neural mechanisms underlying a broad range of natural behaviors in small animals.

Hybridoma screening in semi-solid medium typically employs antigens labeled with visible fluorophores (e.g., FITC, AF488) to enable single-step identification of antibody-secreting clones. However, conventional chemical conjugation via NHS-esters or isothiocyanate groups frequently modifies lysine residues located within epitopes, potentially abrogating antibody recognition of these critical regions. Here, we describe a SpyTag SpyCatcher-based site-specific labeling strategy that circumvents epitope damage during semi-solid medium screening. A 16-amino-acid SpyTag was genetically fused to the C-terminus of the target antigen, enabling covalent conjugation to an sfGFP SpyCatcher fluorescent probe. In semi-solid medium supplemented with SpyTag-antigen and sfGFPSpyCatcher, positive hybridoma clones were readily identified by distinct fluorescent halos, whereas negative clones showed no detectable signal. Notably, the site-specific method yielded a significantly higher frequency of fluorescence-positive clones compared to the conventional AF488-labeled antigen method, suggesting that epitope preservation enhances screening recovery. Furthermore, this approach did not impair hybridoma growth or final clone positivity, offering a simple, rapid, and epitope-compatible method for monoclonal antibody screening.

Molecular staging of Alzheimer's disease (AD) increasingly defines transition boundaries along single-cell pseudo-progression trajectories, yet whether such boundaries reproduce across brain regions, cohorts and molecular modalities is rarely tested. We present a permutation-controlled audit that combines nine boundary-detection algorithms with a fixed marker panel and four orthogonal reproducibility axes-algorithmic consensus, region, cohort and modality. On synthetic data with planted ground-truth boundaries the audit reaches 100% sensitivity and 94% specificity, rejecting four distinct artefact classes each by a different axis. Applied to the Seattle Alzheimer's Disease Brain Cell Atlas middle temporal gyrus, it localizes a transition that is robust across algorithms and recovered in most cell types but does not generalize: its leading marker is attenuated or absent in prefrontal cortex, entorhinal cortex and cerebrospinal fluid, and an apparent cross-region conservation of glial metabolic genes proves to be a global-expression offset rather than a shared program. The same audit nonetheless certifies an externally validated marker (astrocytic PTGDS) as reproducible across regions and modalities, showing that it separates generalizable anchors from dataset-specific ones rather than rejecting all signals. We provide this four-axis audit as a transferable, code-available standard to apply before a trajectory boundary is read as a biological stage, in AD and other progressive proteinopathies.

Linker chemistry and conformation are central determinants of PROTAC activity, shaping ternary-complex geometry, cooperativity, target-lysine presentation and cellular permeability. Existing linker generators often lack explicit control over linker flexibility, require predefined attachment sites and linker lengths, or produce structures that demand substantial geometric correction, limiting their utility in practical PROTAC design. Here we introduce FlexiTAC, a Bayesian flow network that jointly generates linker atom types and coordinates from the warhead and E3-ligase-ligand contexts. We also assemble PROTAC-3D, a quality-controlled collection of 63,554 component-resolved PROTAC structures for model training, and PROTAC-Bench, which covers molecular quality, fragment preservation, geometric fidelity, conformational stability, fragment awareness, rediscovery and sampling efficiency. Compared to the best 3D baseline models, FlexiTAC improves validity by 12.0-12.7% and achieves the highest PoseBusters pass rate of 79.5%-80.0%. A differentiable guidance module shifted generated linkers along a conformational ensemble-derived rigidity axis without retraining the generator. In silico case studies further show that the model can accept crystal-derived, redocked or predicted structural inputs. Together, FlexiTAC, PROTAC-3D and PROTAC-Bench establish an integrated and reproducible framework for data-driven PROTAC linker design, combining controllable structure-conditioned generation with standardized training data and evaluation protocols. This framework expands the linker chemical and conformational space accessible to computational exploration, provides a foundation for future method development and enables the systematic generation of structure-conditioned linker designs with tunable conformational flexibility.

Abstract Imaging-based spatially resolved transcriptomics (imSRT) technologies provide high-throughput molecular-resolution spatial characterization of genes within cells. Conventional analysis methods to identify cell-types and states in imSRT data rely on gene count matrices derived from tallying the number of mRNA molecules detected for each gene per segmented cell, thereby overlooking subcellular heterogeneity that can be useful in defining cell states. To take advantage of the molecular-resolution information in imSRT data and potentially identify cell-states based on subcellular heterogeneity, we developed STARIT (Spatial Transcriptomics As Rasterized Image Tensors). STARIT converts transcripts within segmented cells in imSRT data into an image-based tensor representation that can be combined with deep learning computer vision models for downstream analysis. Using simulated and real imSRT data, we demonstrate that STARIT distinguishes transcriptionally distinct cell-types and further separates cell states based on subcellular transcript localization, which conventional gene count analysis fails to capture. By providing a standardized framework to encode subcellular molecular information in imSRT data, STARIT will enable deeper insights into subcellular heterogeneity and enhance the identification and characterization of cell-types and states that are overlooked by gene count representations.

Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis have shown great promise in treating bladder cancer and are now part of the standard treatment for advanced disease. However, many patients still fail to respond to treatment and at present many biomarkers are assessed but have yet shown only limited results. Therefore, with the advent of combination treatments and the increase of immune related adverse event, the search for reliable predictive biomarkers is paramount. Using a multiplexed enhanced FRET-FLIM based technique (QF-Pro) we quantified the interaction of PD-1/PD-L1, CTLA-4/CD80 and TIGIT/CD155 immune checkpoints in a pre-treatment TMA of 46 patients treated with atezolizumab. The association between higher PD-1/PD-L1 ICP interaction state and treatment efficacy was demonstrated in the male sample cohort, where it identified patients with better PFS. Conversely, patients exhibiting higher CTLA-4/CD80 engagement had a worse response to atezolizumab. Remarkably, the dual assessment of patients with high PD-1/PD-L1 and low CTLA-4/CD80 allowed to identify the best responders. These results indicate that the monitoring of patients immune profile in urothelial carcinoma might be critical in identifying patients who may benefit from combination therapy.

Proteins are dynamic molecules existing in diverse conformational states underlying their biological functions. Although recent approaches have enabled diverse conformational sampling by emulating molecular dynamics simulations, perturbing evolutionary information, or steering internal mechanisms of structure prediction models, predicting conformations resulting from major domain motions or motions that occur over long timescales still remains a challenge. To this end, we introduce GNMCADS, a conformational sampling strategy that enhances the diversity of protein diffusion models by selectively annealing the conditioning signal guided by the intrinsic dynamical organization of the sampled protein. Further, we implement GNMCADS in the diffusion module of AlphaFold3, enabling the generation of diverse protein conformations. When benchmarked across 92 proteins that include 54 class A GPCRs, 15 transporters, and 23 proteins with major domain movements, GNMCADS exhibits improved sampling diversity compared to other current conformational sampling methods.

Methods for installing synthetic functions on living cell surfaces provide powerful approaches for imaging, sensing, and manipulating cell behavior, but many require genetic modification of the target cell or chemical modification of the plasma membrane. Here, we repurpose the glycosylphosphatidylinositol-anchored protein (GPI-AP)-binding toxin aerolysin as a modular chassis for non-genetic cell-surface functionalization. We show that a non-cytotoxic, monomeric aerolysin mutant retains high-affinity and GPI-AP-dependent cell binding when genetically fused to diverse protein cargos. Fluorescent protein-aerolysin fusions robustly label multiple cell types and remain predominantly associated with the cell surface for at least 24 h, in contrast to wheat germ agglutinin, which is extensively internalized. Aerolysin can also be equipped with SpyTag/SpyCatcher to enable modular assembly with independently expressed protein cargos. Importantly, aerolysin supports functional rather than solely optical modification of the cell surface: fusion to the proximity-labeling enzyme APEX2 enables extracellular protein biotinylation, while fusion to HUH endonuclease tags enables covalent attachment of synthetic DNA to living cells. Using this latter architecture, we developed a DNA hairpin sensor that converts cell-surface nuclease activity into a fluorescent signal and distinguishes cells with different levels of extracellular nuclease activity. Together, these results establish non-cytotoxic aerolysin as a genetically encoded, soluble adapter for installing proteins, enzymes, and programmable nucleic acids onto living cells without modification of the target-cell genome.

Label-free Raman imaging of living cells is photon limited: at exposures compatible with cellular dynamics, single-pixel spectra carry about one count per channel on a dominant smooth background. We present an unmixing framework in which the decoder of a physics-constrained autoencoder is restricted to a data-driven spectroscopic dictionary: band centers,widths, and pseudo-Voigt shapes are measured from the dataset and fixed, and the network learns only nonnegative band amplitudes, a smooth B-spline background, and a per-pixel gain.First, on slit-scanning images of HeLa cells (532 nm) the dictionary yields spike-free component spectra that read as band tables, including a resonance-enhanced cytochrome-c-associated component matching literature spectra, and the most stable decomposition against the component number. Second, the dictionary and initialization calibrated at 1 s exposure perline transfer to 100 ms per line (12 s sweeps): cytochrome-c spectral identity survives a single sweep (correlation 0.92) while its map remains photon limited; the dictionary provides spectral physicality, and the transferred initialization prevents a structural collapse that global map correlations miss; in a measurement-derived phantom the dictionary estimator holds thecytochrome-c spectrum to 17-19{degrees} spectral angle at 100 ms, where classical factorizations and free decoders lose it (55-64{degrees}). Estimation on the count-equivalent detector output uses a calibrated shifted-Poisson quasi-likelihood. Third, evaluation must be time matched:correlation against a separately acquired reference saturates through slow specimen drift and acquisition mismatch rather than photon noise, and the self-consistency of learned denoisers is inflated by shared bias; time-matched self-consistency and independent cross-checks areproposed.

Constructing microbial association networks is a common strategy for exploring relationships among taxa in microbiome studies. Although marginal correlation methods are easy to implement and allow formal inference, they can produce spurious edges driven by indirect associations through other taxa. Conditional graphical-modeling methods aim to recover direct associations, but many rely on Gaussian or linear assumptions and often provide limited uncertainty quantification. We propose a conditional, nonparametric approach based on the scaled expected conditional covariance (SEcov). SEcov measures population-level conditional association by residualizing each taxon with respect to the remaining taxa and scaling the resulting expected conditional covariance. The resulting estimator can incorporate flexible machine-learning methods for conditional-mean estimation and admits asymptotic normal inference, enabling p-values and confidence intervals for taxon-pair associations. We demonstrate through simulation studies that our proposed approach improves network recovery relative to other methods, and we illustrate the new method via construction of a co-occurrence network for the vaginal microbiome during pregnancy. IMPORTANCEHigh-throughput sequencing has made it possible to characterize microbial communities at large scale, and network analysis is widely used to summarize relationships among taxa. However, networks based on marginal correlations may include indirect associations, whereas many conditional graphical models rely on assumptions that may be difficult to justify for sparse, zero-inflated, compositional microbiome data. SEcov offers a practical alternative by estimating conditional associations nonparametrically and attaching inferential uncertainty to individual edges. This allows investigators to construct microbiome networks using statistically interpretable evidence for taxon-pair associations, rather than relying solely on arbitrary correlation cutoffs or regularization tuning parameters.

Lysine acylations, including lactylation (Klac), are pivotal regulators of cellular physiology. However, their analysis is currently bottlenecked by antibody enrichment strategies that suffer from sequence bias and require milligram-scale protein inputs, severely precluding the profiling of scarce clinical biopsies and rare cell populations. Here we present ChemIntelligence, an acyl-NHS chemistry-empowered derivatization strategy that rapidly generates unprecedented acylation-specific spectral libraries, exemplified by over 2.5x10^9 human Klac peptides, enabling cross-species reference atlases. Integrated with Prosit-based rescoring, these libraries substantially increase Klac identifications across diverse proteomic datasets. Leveraging this spectral resource, we devised ChemIntelligence Scope, a reproducible, multiplexed parallel reaction monitoring (PRM) platform that quantifies hundreds of Klac peptides per injection from as little as ~200 ng of cell lysates, clinical biopsies, and even true single cells - revealing functional Klac signatures inaccessible to conventional methods. The ChemIntelligence pipeline also extends seamlessly to lysine nicotinylation, underscoring its broad adaptability for discovering and profiling new acylations. Together, these chemical and computational advances establish a scalable, antibody-free framework for acyl-proteome mapping that overcomes input constraints and enables deep functional insights from otherwise intractable biological samples.

Spatial-omics workflows enable molecular analysis within tissue spatial context. Despite the prognostic value of tissue stiffness, these approaches have not incorporated direct, mechanical measurements. This omission reflects several challenges, including sample requirements, low throughput, specialized equipment, and complex data registration. Here, we introduce mechanoMaST (mechanics mapped to spatial transcriptomics), the first workflow to combine absolute mechanical measurements with spatial-omics. It pairs atomic force microscopy-based nanoindentation stiffness maps with spatial transcriptomics maps from adjacent tissue cryosections. The two modalities are then computationally co-registered to enable direct spatial correlation at 100 um resolution, with mapping accuracy quantified through error propagation, providing ground-truth mechanical data directly linked to spatial gene expression. We demonstrate mechanoMaST in human colorectal cancer liver metastasis, generating a spatial resource from 10 patients and revealing a four-gene stiffness signature. mechanoMaST is readily adaptable to other tissues across development and disease, and extendable to additional spatial-omics modalities in adjacent sections.

Genome-wide coexpression analysis of intact tissue samples is a powerful approach for identifying reproducible signatures of cell types and states, since it can survey vast numbers of individuals, cells, and transcripts. However, it can be difficult to optimize gene coexpression network construction and compare results from independent analyses. To address these challenges, we developed OMICON (theomicon.ucsf.edu) for research on human brain gene coexpression networks. OMICON contains gene expression data from >17K normal and neoplastic human brain samples with standardized metadata. Systematic analysis of independent datasets identified >250K gene coexpression modules, which were characterized and compared via enrichment analysis with >40K gene sets. All modules are discoverable via an advanced search engine that can filter by genes, metadata, and enrichment results. Analyses can also be browsed with an interactive workflow visualization tool, and users can communicate within OMICON using @mention functionality to support communal research on human brain gene coexpression networks.